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Abstract : We consider the two-dimensional problem of recovering globally in time the heat distribution 
on the surface of a layer inside of a heat conducting body from two interior temperature measurements. The 
problem is ill-poscd. The approximation function is represented by a two-dimensional Sine series and the 
error estimate is given. 
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1. Introduction 

In this paper, we consider the problem of recovering the heat distribution on the surface of a thin layer 
inside of a heat conducting body from transient temperature measurements. The problem is raised in many 
applications in Physics and Geology. In fact, in many physical situation (see, e.g. [B]) we cannot attach 
a temperature sensor at the surface of the body (for example, the skin of a missile). On the other hand, 
we can easily measure the temperature history at an interior point of the body. Hence, to get the heating 
history in the body, we want to use temperature measured in the interior of the body. In the simplest model, 
the heat-conducting body is assumed to have a constant conductivity and represented by the half-line a; > 
(see, e.g. [C, EM, LN, TV]),. While giving many useful results, this model is not suitable for the case of a 
body having a series of superposed layers, each of which has a constant conductivity. 

Precisely, we shall consider the problem corresponding to a thin layer of the body represented by the strip 
Rx (0, 2), say. Let u be the temperature in the strip. For the uniqueness of solution, we shall have to measure 
the temperature history at two interior lines R x {y = 1} and R x {y = 2}. From these measurements, we 
can identify uniquely the heating history inside of the layer (see, e.g., [B]). The problem is of finding the 
surface heat distribution u{x,0,t) = v{x,t). In fact, despite uniqueness, the global solution in L^(R x R_|_) 
is unstability and hence, in this point of view, a sort of regularization is in order. 

As discussed in the latter paragraph, the main purpose of our paper is to present a regularization of the 
problem. Moreover, an effective way of approximating the heat distribution v is also worthy of considering. 
There are many methods for regularizing the equation (see [TA, B]). In the most common scheme (see, e.g., 
[Blackwell] ) , the computation is divided into two steps. In the first step, one considers the problem of finding 
the heat flux Uy{x, 1, t) from the interior measurements u{x, 1, t), u{x, 2, t). The problem is classical and can 
be changed to the one of finding the solution of a convolution of Volterra type which can be solved in any 
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finite time interval by tiie iteration (see, e.g., [F]). But even in the "classical" problem as mentioned in many 
documents, it is worthy of insisting that the problem is ill-posed if we consider the problem over the whole 
time interval R+ with respect to the L^-norm and the literature on this way is very scarce. In the second 
step, one considers the "really" ill- posed problem which is of recovering the surface temperature history 
u{x,0,t) from data u{x,l,t),Uy{x,l,t). In the present paper, we can (and shall) regularize the function 
u{x,0,t) directly from u{x,l,t),u{x,2,t) without using the flux function Uy{x,l,t). We emphasize that, 
using our method, we can unify two steps of the classical scheme and find simultaneously two functions 
u{x,0,t) and Uy{x,l,t). However, our main purpose is of regularizing the surface temperature. Hence, we 
shall omit the problem of finding the interior heat flux Uy{x, 1, t). 

Moreover, using the method of truncated integration, one can approximate the Fourier transform of the 
solution by a function having a compact support in R^. Therefore, the solution can be represented by an 
expansion of two dimensional Sine series (see [AGTV]). The Sine method is based on the Cardinal functions 

^/ ,^/ N sin\(Tr(z — pd))/d] „ , 

^(^'^)(-)= .(.-Jvi 

which dates back to the works of many mathematicians (Bohr, de la Vallee Poussin, E. T. Whittaker, 
...). The one dimensional version of the method is studied very clearly and completely in [S]. Some primary 
results related to the two dimensional Sine approximation were given in [AGTV]. As is known, the Sine series 

1/2 

converges very rapidly at an incredible 0(e~'^" ) rate, where c > and n is the dimension of approximation 
(see [S]). Hence, this method, which is new in our knowledge, is very effective. 

The remainder of the present paper is divided into two sections. In Section 2 , we state precisely the 
problem, change it into an integral equation of convolution type, and state main results of our paper. In 
Section 3, we give the proof of main results. 

2. Notations and main results 

Consider the problem of determining the heat distribution 

u{x,0,t) = v{x,t) (1) 

where u satisfies 

du 

Au - — = a; e R, < y < 2, t > 0, (2) 
subject to the boundary conditions 

u{x, 2, t) = g{x, t), a: e R, t> 0, (3) 

u{x,l,t) = f{x,t), X e R, t>0, (4) 



and the initial condition 



u{x, y, 0) = 0, a; e R, < y < 2. (5) 



Here f,g are given. We shall transform the problem (l)-(5) into a convolution intergral equation. 
Put 

r{x, y, t, r) = exp ^ j (6) 
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and 



We have 



G{x, y, t, ^, T], t) = T{x, y, t, ^, rj, r) - T{x, 4:-y,t, ^, rj, r). 



+ Grir) + Gt = 0. 



Integrating the identity 

div{GS/u - uVG) - ^{uG) = 

OT 

over the domain R x (1, 2) x {0,t — e) and letting e ^ 0, we have 

+ 00 t +00 t 



j j g{i,T)Gr,{x,y,t,i,2,T)didT + j j Gix,y,t,i,l,T)uy{^,l,T)d^dT 

-oo 

- J J f{i,T)Gn{x,y,t,^,l,T)d^dT + u{x,y,t) 



-oo 



Hence 



+CXD t 



j j G{x, y, t,^,l, T)uy{^, 1, T)d£,dT = -u{x, y, t) + 

-oo 

+ 00 t +00 t 

J J Gr,{x,y,t,^,l,T)fi^,T)d^dT- J J g{^,T)Gr,{x,y,t,^,2,T)d^dT. 



-oo 



Letting ?/ — > 1+ in (8), we have 

+ 00 t 



1 f {x-0' 

exp ' 



-oo 



2TT{t - t) 



-oo 



1 / (2^-0' + 4 

exp ' 



+00 t 

-fix,t)~— / / - — 3^ exp 



A{t -t) J 27r(t - r) V 4(f - r) 

(x -0^ + 4 



27r y J {t- t)2 

-00 



4(i-T) 



fi^,r)d^dT 



+ 



+ C50 t 

ill 



9{^,r)-^e^p (J-^-Jl±l) d^dr. 



-00 



4{t-T) 



We put N{x, y, t, r;, r) = T{x, y, t, S,, r], r) - r(x, -y, t, ^, r), r) 
Integrating the identity 

d 

div{NVu - uVN) - ^("A'') = 



over the domain (— n, n) x (0, 1) x (0, i — e) and letting n — ^ oo, e — ^ 



+00 t 



+ 00 t 



j j N{x,y,t,^,l,T)Uy{^,l,T)d^dT - j j f{^,T)Nr,{x,y,t,^,l,T)d^dT 

-oo 

T)Nn{x, y, t, ^, 0, T)d^dT - u{x, y, t) = 0. 



-oo 

-t-oo t 



-oo 



Letting ?/ — > 1 , the identity (10) becomes 



1 



+OC' t 



1 



27r y J t-T 

-oo 

-t-oo t 



exp 



4(i-r) 



exp 



(x-g)^+4 
4(t-r) 



1 
"2^ 



-oo 
-t-oo t 



+ 



t)2 '^^^ ' 



(x -0^+4 

4(t-T) 



d^dr 



^//(^^"P(" 

-oo 



(a;-0' + l 



4(i-T) 



v{^,T)d^dT -3f{x,t) =0 



From (9) and (11), we have the main convolution equation 



TT J J {t-Tf 

-oo 



exp - 



(a;-0'+4 



/(C,T)dedr 



4(i-T) 

+ ^//^(^-)(*^-p(-^%^)^^- 

(x - 0' + 1 



+00 i 



-oo 

+00 t 



2tt J J {t- r)2 

-oo 



exp 



4(i-T) 



v{^,T)d^dT-4f{x,t) = 



which can be rewritten as 

S * v{x, t) =2R* f{x, t)-S* g{x, t) + 4/(x, t) 
where we define that v{x,t) = f{x, t) — g{x,t) = as t < 0, 



R{x,t) 



F (- 



At 







{x,t) e R X [0, +oo) 
{x,t) e R X (-00,0) 



and 



S{x,t) 



^f^exp(-4ti) (x,i)eRx [0,-Koo) 
\ (a;,t) e R x (-00,0) 



Put F{x, t) = 2R* f{x, t)- S* g{x, t) + Af{x, t). 
Taking the Fourier-transform of (12), we get 

S{z, r)v{z, r) = F{z, r) 



where 



+00 +00 



S{z,r) 



= ^ J J S{z,r)e-''-'''+*''Uxdt 



2e 



COS -j= \/ y/ z"^ + — — isgn{r) sin \/ \/ z^ + 



and 



S{z,r) 



2e 
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We have 
Theorem 1 

Let 7 G (0,2) and e € (o,e^). 

Assume that Vq G L^(R^) is the (unique) solution of (12) corresponding to the exact data fo, go S i^(R^) 
and that f,g G L^(R'^) are measured data satisfying \\f — /0II2 < e and — go\\2 < e wiere ||.||2 is tiie 
i2(R^)-norm. 

Then we can construct from g,f a function G L"^ such that 

< VCs^-'' + vie) 

where C is constant and rj(e) as e J, 0. 

Moreover, if we assume in addition that vq G -ff'"(R^) fl L^(R^), m > and 
< e < min{e~^ , e~^™ } then 



-wolla < -D ( In ( - 



where D > depends on vq. 
Theorem 2 

With Ve as in theorem 1, we have 



wiiere 



Vs{x,t)= ^ vJ—,— \S{m,n/as){x)S{n,n/ae){t) 

n=-oo|„|<|„| 



„, , , sin\-K{z — vd\ I (S\ „ , 



■k{z — pd)/d 



3. Proofs 

Proof of theorem 1 

We put 



and 



F{x, t) = 2R* fix, t)-S* g{x, t) + 4f{x, t) 
Fo{x, t)=2R* fo{x, t)-S* goix, t) + 4fo{x, t) 
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then 



Put 



< {4 + 2\\R\\,)\\f-foh + \\S\\,\\g-goh 

< {A + 2\\R\\, + \\S\\,)e. 



27r J S{z, r) 



where De = {{z,r)/ \z\ < be and \r\ < b^} with be = ^^^^^ In ^• 



We have 



||Ve-t^0||2 = ll'>£-^'0||2 



-I 

< £2-7 



F{z,r)-Fo{z,r) 



dzdr + 



dzdr 



S{z,r) 

{A + 2\\R\\, + \\S\\{f+ j \vo{z,r)\' dzdr 



/ 



If we put r7(e) = / \vo{z, r)!^ dzdr then r]{e) ^ as e | 0. 

Now, we assume vq e H'^{B?), m > 0, put 

x/2 



In ( . li^, . I > 1 



and 



Vln'"(l/e 



Ve{x,t) = — I l^^e'^^'+'^^dzdr 
J S{z,r) 



(13) 



We have 



||t^e-t^o||2 = j 



F{z, 


r)-Fo{z,r) ' 




S{z,r) 


2 



-dzdr + j \vo{z,r)\^ dzdr 



(2^ + H)"|£.o(z,r) 

(z2 + 7-2)™ 



where 



< 4£2(4 + 2||ii||i + ||5||i)2e^^^«=+ y 

R'\Q. 



Ci=max|4(4 + 2||ii||i + ||5||i)\||(^2^r2r/%(^,r)|[2'"| 



-dzdr 
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This implies that 



Ibe-volla < Ci 
< C, 



l/e 



ln2"(l/£) ■ ln""(l/£) 



2m/ 



ln2"(l/e) 



where D = ^Ci(l + 2™). 
This completes the proof. 
Proof of theorem 2 
We have 



supp Ve C Ds c [-ae,ae] X [-ae,as]. 



As in [AGLT], p.l21, we have theorem 2. 
This completes the proof. 
4. Numerical results 

We present some results of numerical comparison of the regularized representation of the solution given 
by theorem 2 and the corresponding exact solution of the problem. 

Let the problem 



Am- — =0, e Rx (0,2), t >0 



1 -rr^-l 1 

u{x,l,t) = je ** ; u{x,2,t) = -e •,u{x,y,0) = 



whose the unknown is 



(14) 
(15) 

(16) 



v{x,t) = u{x,0,t), X > xq > 0, t>0 
The exact solution of this problem is 

1 ^ 
v{x,t) = -e « . 

The approximated solution is calculated from the expansion of two-dimensional Sine series given by theorem 
2 associated to formula (13) in which 



F = 4- 



S = 2e-^'^ 



Thus we have 



F-Fr, 



z,2(R') 



= £ 
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which is a perturbation similar to the one operated on the two given functions / and g. 

With £ = iV = 50 (the size of the double series) and for {x, t) e [0.25, 1.3] x [0, 4] we have drawn the 
corresponding approximate surface solution {x,t) — > Ve{x,t) in Fig.l. 

To calculate the double integral in (10) we have used the rectangle rule which gives good accuracy if 
one integrates on the interval [e, 1/e] x [e, 1/s]. The time of calculus with a good computer is very long: 2 
hours for 900 points M — {x,t) e [0.25, 1.3] x [0,4]. It is the reason for which we are limited ourselves to a 
relatively small size of the double series {N = 20). For comparison in Fig. 2 we have drawn the exact solution 
{x, t) — > u(x, t). 

approximate solution 




Fig.l: regularized solution of the problem (14), (15) 



exact soiution 




Fig.l: exact solution of the problem (14), (15) 

Using the same method as previously we have drawn in Fig. 3 the surface {x, t) E [0, 1] x [0, 4] 
which is the regularization of the following problem 



Ve{x,t) 



Au-^=0, (x, y) e R X (0, 2),t>0 



(17) 
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u{x,l,t)=0; u{x,2,t) = -e—n—; u{x,y,0) = (18) 
the unknown being v{x, t) = u{x, 0, t). The exact sohition 




is represented in Fig. 4. 



approximate solution 




Fig. 3: regularized solution of the problem (17), (18) 

exaot solution 




Fig. 4: exact solution of the problem (17), (18) 
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